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Heavy Ion Collisions (HIC) represent a unique tool to probe the in-medium nuclear 
interaction in regions away from saturation and at high nucleon momenta. In this report 
we present a selection of reaction observables particularly sensitive to the isovector part 
of the interaction, i.e. to the symmetry term of the nuclear Equation of State (EoS) At 
low energies the behavior of the symmetry energy around saturation influences dissipation 
and fragment production mechanisms. Predictions are shown for deep-inelastic and frag- 
mentation collisions induced by neutron rich projectiles. Differential flow measurements 
will also shed lights on the controversial neutron/proton effective mass splitting in asym- 
metric matter. The high density symmetry term can be derived from isospin effects on 
heavy ion reactions at relativistic energies (few AGeV range), that can even allow a "di- 
rect" study of the covariant structure of the isovector interaction in the hadron medium. 
Rather sensitive observables are proposed from collective flows and from pion/kaon pro- 
duction. The possibility of the transition to a mixed hadron-quark phase, at high baryon 
and isospin density, is finally suggested. Some signatures could come from an expected 
"neutron trapping" effect. 



1. Introduction 

The symmetry energy E sym appears in the energy density e(p, p 3 ) = e(p)+pE sym (p 3 / p) 2 + 
0(p3,/ p) 4 + .., expressed in terms of total (p = p p + p n ) and isospin (p 3 = p p — p n ) densities. 
The symmetry term gets a kinetic contribution directly from basic Pauli correlations and 
a potential part from the highly controversial isospin dependence of the effective interac- 
tions [ 1J. Both at sub-saturation and supra-saturation densities, predictions based of the 
existing many-body techniques diverge rather widely, see [ |2] . We take advantage of new 
opportunities in theory (development of rather reliable microscopic transport codes for 
HIC) and in experiments (availability of very asymmetric radioactive beams, improved 
possibility of measuring event-by-event correlations) to present results that are severely 
constraining the existing effective interaction models. We will discuss dissipative colli- 
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Figure 1. Distribution of the octupole moment of primary fragments for the 132 Sn+ Ni 
reaction at 10 AMeV (impact parameters (a):6 = 6/m, (b):7/m, (c):8/m). Solid lines: 
asysoft. Dashed lines: asystiff 

sions in a wide range of energies, from just above the Coulomb barrier up to a few AGeV. 
The transport codes are based on mean field theories, with correlations included via hard 
nucleon-nucleon elastic and inelastic collisions and via stochastic forces, selfconsistently 
evaluated from the mean phase-space trajectory, see [HJ ESSIE]. Stochasticity is essential 
in order to get distributions as well as to allow the growth of dynamical instabilities. 

2. Isospin effects on Deep-Inelastic Collisions 

Dissipative semi-peripheral collisions at low energies, including binary and three-body 
breakings, offer a good opportunity to study phenomena occurring in nuclear matter under 
extreme conditions with respect to shape, excitation energy, spin and N/Z ratio (isospin). 
In some cases, due to a combined Coulomb and angular momentum (deformation) effect, 
some instabilities can show up This can lead to 3-body breakings, where a light cluster 
is emitted from the neck region. Three body processes in collisions with exotic beams 
will allow to investigate how the development of surface (neck-like) instabilities, that 
would help ternary breakings, is sensitive to the structure of the symmetry term around 
(below) saturation. In order to suggest proposals for the new RIB facility Spiral 2, [ 
[7] we have studied the reaction 132 Sn + 64 Ni at lOAMeV in semicentral events, impact 
parameters b = 6,7, 8/m, where one observes mostly binary exit channels, but still in 
presence of large dissipation. Two different behaviors of the symmetry energy below 
saturation have been tested: one (asysoft) where it is a smooth decreasing function 
towards low densities, and another one (asystif f) where we have a rapid decrease, [ 
P. The Wilczynski plots, kinetic energy loss vs. deflection angle, show slightly more 
dissipative events in the asystiff case, consistent with the point that in the interaction 
at lower densities in very neutron-rich matter (the neck region) we have a less repulsive 
symmetry term. In fact the neck dynamics is rather different in the two cases, as it can be 
well evidenced looking at the deformation of the PLF/TLF residues. The distribution 
of the octupole moment over the considered ensemble of events is shown in FigC] for 
the three considered impact parameters. Except for the most peripheral events, larger 
deformations, strongly suggesting a final 3-body outcome, are seen in the asystiff case. 
Now, due to the lower value of the symmetry enrgy, the neutron-rich neck connecting 
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the two systems survives a longer time leading to very deformed primary fragments, 
from which eventually small clusters will be dynamically emitted. Finally we expect 
to see effects of the different interaction times on the charge equilibration mechanism, 
probed starting from entrance channels with large N/Z asymmetries, like 132 Sn(N/Z = 
1.64) + 58 Ni(N/Z = 1.07). Moreover the equilibration mechanism is also directly driven 
by the strenght of the symmetry term. 

3. Isospin Dynamics in Neck Fragmentation at Fermi Energies 

It is now quite well established that the largest part of the reaction cross section for 
dissipative collisions at Fermi energies goes through the Neck Fragmentation channel, with 
IMFs directly produced in the interacting zone in semiperipheral collisions on very short 
time scales [|H]. We can predict interesting isospin transport effects for this new frag- 
mentation mechanism since clusters are formed still in a dilute asymmetric matter but 
always in contact with the regions of the projectile-like and target-like remnants almost at 
normal densities. Since the difference between local neutron-proton chemical potentials 
is given by \i n — fi p = ^E sym (p 3 / p) , we expect a larger neutron flow to the neck clusters 
for a stiffer symmetry energy around saturation, [ |U] . The isospin dynamics can be 
directly extracted from correlations between N/Z, alignement and emission times of the 
IMFs. The alignment between PLF — IMF and PLF — TLF directions represents a 
very convincing evidence of the dynamical origin of the mid-rapidity fragments produced 
on short time scales [ E] . The form of the <& p i ane distributions (centroid and width) 
can give a direct information on the fragmentation mechanism [E]- Recent calculations 
confirm that the light fragments are emitted first, a general feature expected for that 
rupture mechanism [E|- The same conclusion can be derived from direct emission time 
measurements based on deviations from Viola systematics observed in event-by-event ve- 
locity correlations between IMFs and the PLF /TLF residues [El El El- We can figure 
out a continuous transition from fast produced fragments via neck instabilities to clusters 
formed in a dynamical fission of the projectile(target) residues up to the evaporated ones 
(statistical fission). Along this line it would be even possible to disentangle the effects of 
volume and shape instabilities. A neutron enrichment of the overlap ("neck") region is 
expected, due to the neutron migration from higher (spectator) to lower (neck) density 
regions, directly related to the slope of the symmetry energy [E]- A very nice new anal- 
ysis has been presented on the Sn + Ni data at 35 AMeV by the Chimera Collab., Fig. 2 
of ref. [Q3j. A strong correlation between neutron enrichemnt and alignement (when the 
short emission time selection is enforced) is seen, that can be reproduced only with a stiff 
behavior of the symmetry energy. This is the first clear evidence in favor of a relatively 
large slope (symmetry pressure) around saturation. 

4. Effective Mass Splitting and Collective Flows 

The problem of Momentum Dependence in the Isovector channel (Iso — MD) is still 
very controversial and it would be extremely important to get more definite experimen- 
tal information, see the recent refs. [ El El El El El E] • Intermediate energies are 
important in order to have high momentum particles and to test regions of high baryon 
(isoscalar) and isospin (isovector) density during the reactions dynamics. Collective flows 
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[ |2J] are very good candidates since they are expected to be very sensitive to the mo- 
mentum dependence of the mean field, see [|22J[T]. The transverse flow, Vi(y,p t ) = 
provides information on the anisotropy of nucleon emission on the reaction plane. Very 

important for the reaction dynamics is the elliptic flow, V^iy^Pt) — ( Px /" )• The sign 
of V2 indicates the azimuthal anisotropy of emission: on the reaction plane (V2 > 0) or 
out-of-plane (squeeze — out, V 2 < 0) [ I2T] 122] . We have then tested the I so — MD of 
the fields just evaluating the Difference of neutron/proton transverse and elliptic flows 
^12 \ViPt) = Vi 2 (y,Pt) — Vi^iUiPt) at various rapidities and transverse momenta in 
semicentral (b/b max = 0.5) 197 Au+ 197 Au collisons at 250AMeV, where some proton data 
are existing from the FOPI collaboration at GSI [EH] IS]. The transport code has been 
implemented with a BGBD — like [ I2J3 12^] mean field with a different (n, p) momentum 
dependence, see [[THJ [T7J HB|, that allow to follow the dynamical effect of opposite n/p ef- 
fective mass splitting while keeping the same density dependence of the symmetry energy. 
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Figure 2. Difference between proton and neutron V\ flows in a semi-central reaction 
Au+Au at 250 AMeV for three rapidity ranges. Upper Left Panel: \y^\ < 0.3; Upper 
Right: 0.3 < \y^\ < 0.7; Lower Left: 0.6 < \y^\ < 0.9. Lower Right Panel: Comparison 
of the V\ proton flow with FOPI data [|2S1 for three rapidity ranges. Top: 0.5 < < 
0.7; center: 0.7 < \y^\ < 0.9; bottom: 0.9 < \y^\ < 1.1. 

For the difference of nucleon transverse flows, see Fig. |2J the mass splitting effect is 
evident at all rapidities, and nicely increasing at larger rapidities and transverse mo- 
menta, with more neutron flow when m* < m*. Just to show that our simulations give 
realistic results we compare in lower right panel of Fig. |2J with the proton data of the 
FOPI collaboration for similar selections of impact parameters rapidities and transverse 
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momenta. The same analysis has been performed for the difference of elliptic flows, [ 
IT7] . Again the mass splitting effects are more evident for higher rapidity and tranverse 
momentum selections. In particular the differential elliptic flow becomes negative when 
m* < m*, revealing a faster neutron emission and so more neutron squeeze out (more 
spectator shadowing). The measurement of n/p flow differences appears essential. Due to 
the difficulties in measuring neutrons, our suggestion is to measure the difference between 
light isobar flows, like 3 if vs. 3 He and so on. We expect to clearly see the effective mass 
splitting effects, maybe even enhanced due to larger overall flows shown by clusters, see [ 

HI 123 



5. Relativistic Collisions 

Finally we focus our attention on relativistic heavy ion collisions, that provide a unique 
terrestrial opportunity to probe the in-medium nuclear interaction at high densities. An 
effective Lagrangian approach to the hadron interacting system is extended to the isospin 
degree of freedom: within the same frame equilibrium properties (EoS, [ 28]) and trans- 
port dynamics [ EHl E01 can be consistently derived. Within a covariant picture of the 
nuclear mean field, for the description of the symmetry energy at saturation (04 param- 
eter of the Weizsaecker mass formula) (a) only the Lorentz vector p mesonic field, and 
(b) both, the vector p (repulsive) and scalar 5 (attractive) effective fields [ ED E2] can 
be included. In the latter case the competition between scalar and vector fields leads to 
a stiffer symmetry term at high density [ HHJ E] • The presence of the hadronic medium 
leads to effective masses and momenta M* = M + S s , k*^ = k^ — S M , with S s , S M scalar 
and vector self-energies. For asymmetric matter the self-energies are different for protons 
and neutrons, depending on the isovector meson contributions. We will call the corre- 
sponding models as NLp and NLp5, respectively, and just NL the case without isovector 
interactions. 

For the description of heavy ion collisions we solve the covariant transport equation 
of the Boltzmann type [ I2HJ E0| within the Relativistic Landau Vlasov (RLV) method, 
using phase-space Gaussian test particles [ESI, an d applying a Monte-Carlo procedure 
for the hard hadron collisions. The collision term includes elastic and inelastic processes 
involving the production/absorption of the A(1232MeV A ) and N*(1440MeV) resonances 
as well as their decays into pion channels, [EHE5I- A larger repulsive vector contribution 
to the neutron energies is given by the p-coupling. This is rapidly increasing with density 
when the 5 field is included [EH U]- As a consequence we expect a good sensitivity to 
the covariant structure of the isovector fields in nucleon emission and particle production 
data. Moreover the presence of a Lorentz magnetic term in the relativistic transport 
equation [I2H1EII1II] will enhance the dynamical effects of vector fields [ESI- 

Differential flows will be directly affected. In FigEl transverse and elliptic differential 
flows are shown for the 132 Sn+ 124 Sn reaction at 1.5 AGeV (b = Qfm), [EH]- The effect of 
the different structure of the isovector channel is clear. Particularly evident is the splitting 
in the high p t region of the elliptic flow. In the (p + 5) dynamics the high-p t neutrons 
show a much larger squeeze — out. This is fully consistent with an early emission (more 
spectator shadowing) due to the larger p-field in the compression stage. 
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Figure 3. Differential neutron-proton flows for the 132 Sn + 124 Sn reaction at 1.5 AGeV (b = 
6/m) from the two different models for the isovector mean fields. Left: Transverse Flows. Right: 
Elliptic Flows. Full circles and solid line: NLp5. Open circles and dashed line: NLp. 
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Figure 4. Time evolution of the A ±,0,++ resonances and pions ir^ (left), and kaons (K +, ° 
(right) for a central (6 = fm impact parameter) Au+Au collision at 1 AGeV incident energy. 
Transport calculation using the NL, NLp, NLp5 and DDF models for the iso-vector part of the 
nuclear EoS are shown. 

6. Isospin effects on sub-threshold kaon production at intermediate energies 

Kaon production has been proven to be a reliable observable for the high density EoS 
in the isoscalar sector [ETJEH] Here we show that the K 0,+ production (in particular the 
K°/K + yield ratio) can be also used to probe the isovector part of the EoS. 

Using our RMF transport approach we analyze pion and kaon production in central 
197 Au + 197 Au collisions in the 0.8 — 1.8 AGeV beam energy range, comparing models 
giving the same "soft" EoS for symmetric matter and with different effective field choices 
for E sym Here we also use a Lagrangian with density dependent couplings {DDF, [ 



), recently suggested for better nucleonic properties of neutron stars [EH]- m the DDF 
model the p-coupling is exponentially decreasing with density, resulting in a rather " soft" 
symmetry term at high density. The hadron mean field propagation, which goes beyond 
the "collision cascade" picture, is essential for particle production yields: in particular the 
isospin dependence of the self-energies directly affects the energy balance of the inelastic 
channels. 

Fig. 0] reports the temporal evolution of A ±,0,++ resonances, pions (7r ± '°) and kaons 
(K +, °) for central Au+Au collisions at lAGeV. It is clear that, while the pion yield 
freezes out at times of the order of 50 fm/c, i.e. at the final stage of the reaction (and at 
low densities), kaon production occur within the very early (compression) stage, and the 
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yield saturates at around 20 fm/c. From Fig. 01 we see that the pion results are weakly 
dependent on the isospin part of the nuclear mean field. However, a slight increase 
(decrease) in the 7r~ (7r + ) multiplicity is observed when going from the NL (or DDF) to 
the NLp and then to the NLpS model, i.e. increasing the vector contribution f p in the 
isovector channel. This trend is more pronounced for kaons, see the right panel, due to 
the high density selection of the source and the proximity to the production threshold. 

When isovector fields are included the symmetry potential energy in neutron-rich matter 
is repulsive for neutrons and attractive for protons. In a HIC this leads to a fast, pre- 
equilibrium, emission of neutrons. Such a mean field mechanism, often referred to as 
isospin fractionation [ G] , is responsible for a reduction of the neutron to proton ratio 
during the high density phase, with direct consequences on particle production in inelastic 
iViV collisions. Threshold effects represent a more subtle point. The energy conservation 
in a hadron collision is expressed in terms of the canonical momenta, i.e. for a reaction 
1 + 2 — > 3 + 4 as s in = (k± + k^) 2 = (&3 + k^) 2 = s out . Since hadrons are propagating 
with effective (kinetic) momenta and masses, an equivalent relation should be formulated 
starting from the effective in-medium quantities k*^ = k tJ- — H fl and m* = m+S s , where E s 
and S M are the scalar and vector self-energies. The self-energy contributions will influence 
the particle production at the level of thresholds as well as of the phase space available 
in the final channel. In neutron-rich colliding systems Mean field and threshold effects 
are acting in opposite directions. At low energies, around the production threshold,, the 
energy conservation (i.e. the self energy contributions) is dominant, as we see from Fig. 

in particular for kaons. 

We have to note that in a previous study of kaon production in excited nuclear matter 
the dependence of the K°/K + yield ratio on the effective isovector interaction appears 
much larger (see Fig. 8 of ref.[ 13*1]). The point is that in the non-equilibrium case of 
a heavy ion collision the asymmetry of the source where kaons are produced is in fact 
reduced by the n — > p "transformation" , due to the favored nn — ► pA~ processes. This 
effect is almost absent at equilibrium due to the inverse transitions. Moreover in infinite 
nuclear matter even the fast neutron emission is not present. This result clearly shows 
that chemical equilibrium models can lead to uncorrect results when used for transient 
states of an open system. 

7. Testing Deconfinement at High Isospin Density 

The hadronic matter is expected to undergo a phase transition into a deconfined phase of 
quarks and gluons at large densities and/or high temperatures. On very general grounds, 
the transition critical densities are expected to depend on the isospin of the system, but 
no experimental tests of this dependence have been performed so far. In order to check 
the possibility of observing some precursor signals of some new physics even in collisions 
of stable nuclei at intermediate energies we have performed some event simulations for the 
collision of very heavy, neutron-rich, elements. We have chosen the reaction 238 U + 238 U 
(average proton fraction Z/A = 0.39) at 1 AGeV and semicentral impact parameter 
6 = 7 fm just to increase the neutron excess in the interacting region. After about 
10 fm/c in the overlap region a nice local equilibration is achieved. A rather exotic 
nuclear matter is formed in a transient time of the order of 10 fm/c, with baryon density 



8 



M.Di Toro 



around 3 — 4p , temperature 50 — 60 MeV, energy density ~ 500 MeV fm~ 3 and proton 
fraction between 0.35 and 0.40, likely inside the estimated mixed phase region [ I40j . 

We can study the isospin dependence of the transition densities [EI] m a systematic way. 
Concerning the hadronic phase, we use the relativistic non-linear model of Glendenning- 
Moszkowski (in particular the "soft" GM3 choice) [ I42j . where the isovector part is treated 
just with p meson coupling, and the iso-stiffer NLp5 interaction [EDI- For the quark phase 
we consider the MIT bag model with various bag pressure constants. In particular we 
are interested in those parameter sets which would allow the existence of quark stars [ 
14*3*] . i.e. parameter sets for which the so-called Witten-Bodmer hypothesis is satisfied [ 
El E2]- One of the aim of our work it to show that if quark stars are indeed possible, 
it is then very likely to find signals of the formation of a mixed quark-hadron phase in 
intermediate-energy heavy-ion experiments [ 140] . The structure of the mixed phase is 
obtained by imposing the Gibbs conditions [ES] for chemical potentials and pressure and 
by requiring the conservation of the total baryon and isospin densities 

d H, =d Q \ p<^t,MS) = p<^i>«>), 

Pb = (1-x)Pb +XPb, ft = (l-xK + XPt, (1) 

where x is the fraction of quark matter in the mixed phase. In this way we get the binodal 
surface which gives the phase coexistence region in the (T,p B ,p 3 ) space [HIE EI]. For a 
fixed value of the conserved charge p% we will study the boundaries of the mixed phase 
region in the (T,ps) plane. In the hadronic phase the charge chemical potential is given 
by /i 3 = 2E sym (pB) — ■ Thus, we expect critical densities rather sensitive to the isovector 
channel in the hadronic EoS. 

In Fig. we show the crossing density p cr separating nuclear matter from the quark- 
nucleon mixed phase, as a function of the proton fraction Z/A. We can see the effect of 
the ^-coupling towards an earlier crossing due to the larger symmetry repulsion at high 
baryon densities. In the same figure we report the paths in the (p, Z/A) plane followed 
in the cm. region during the collision of the n-rich 132 Sn+ 132 Sn system, at different 
energies. At 300 AMeV we are just reaching the border of the mixed phase, and we are 
well inside it at 1 AGeV. Statistical fluctuations could help in reducing the density at 
which drops of quark matter form. The reason is that a small bubble can be energetically 
favored if it contains quarks whose Z/A ratio is smaller than the average value of the 
surrounding region [EDI- This corresponds to a neutron trapping effect, supported also by 
a symmetry energy difference in the two phases. In fact while in the hadron phase we have 
a large neutron potential repulsion (in particular in the NLpS case), in the quark phase 
we only have the much smaller kinetic contribution. If in a pure hadronic phase neutrons 
are quickly emitted or "transformed" in protons by inelastic collisions, when the mixed 
phase starts forming, neutrons are kept in the interacting system up to the subsequent 
hadronization in the expansion stage [EHj- Observables related to such neutron "trapping" 
could be an inversion in the trend of the formation of neutron rich fragments and/or of 
the 7r~/7r + , K°/K + yield ratios for reaction products coming from high density regions, 
i.e. with large transverse momenta. 
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Figure 5. Variation of the transition density with proton fraction for various hadronic EoS 
parameterizations. Dotted line: GM3 parametrization; dashed line: NLp parametrization; solid 
line: NLp5 parametrization. For the quark EoS, the MIT bag model with i? 1 / 4 =150 MeV. 
The points represent the path followed in the interaction zone during a semi-central 132 Sn+ 132 Sn 
collision at 1 AGeV (circles) and at 300 AMeV (crosses). 

8. Perspectives 

We have shown that violent collisions of n-rich heavy ions from low to relativistic en- 
ergies can bring new information on the isovector part of the in-medium interaction, 
qualitatively different from equilibrium EoS properties. We have presented quantitative 
results in a wide range of beam energies. At low energies we have shown isospin effects 
on the dissipation in deep inelastic collisions, at Fermi energies the Iso-EoS sensitivity of 
the isospin transport in fragment reactions and finally at intermediate the dependence of 
differential flows on the I so — MD and effective mass splitting. In relativistic collisions we 
have shown the possibility of a direct measure of the Lorentz structure of the isovector 
fields at high baryon density, from differential collective flows and yields of charged pion 
and kaon ratios. Important non-equilibrium effects for particle production are stressed. 
Finally our study supports the possibility of observing precursor signals of the phase 
transition to a mixed hadron-quark matter at high baryon density in the collision, cen- 
tral or semi-central, of neutron-rich heavy ions in the energy range of a few AGeV. As 
signatures we suggest to look at the isospin structure of hadrons produced at high trans- 
verse momentum, as a good indicator of the neutron trapping effect. In conclusion the 
results presented here appear very promising for the possibility of exciting new results 
from dissipative collisions with radioactive beams. 
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